Fractional domain walls from on-site softening in dipolar bosons 
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We study dipolar bosons in a ID optical lattice and identify a region in parameter space — 
strong coupling but relatively weak on-site repulsion — hosting a series of stable charge-density-wave 
(CDW) states whose low-energy excitations, built from "fractional domain walls," have remarkable 
similarities to those of non-Abelian fractional quantum Hall states. Here, a conventional domain 
wall between translated CDW's may be split by inserting strings of degenerate, but inequivalent, 
CDW states. Outside these insulating regions, we find numerous supersolids as well as a superfluid 
regime. The mentioned phases should be accessible experimentally and, in particular, the fractional 
domain walls can be created in the ground state using single-site addressing, i.e., by locally changing 
the chemical potential. 



PACS numbers: 67.85.-d 

I. INTRODUCTION 

The interest in low-dimensional systems, where quan- 
tum fluctuations are particularly strong, has increased 
tremendously during the last decades. Demonstrated 
by theoretical and experimental studies, systems con- 
strained to one or two spatial dimensions exhibit features 
related to statistics and topology which are absent in 
three dimensions. Luttinger liquids, the fractional quan- 
tum Hall effect, anyonic and non-Abelian statistics, and 
graphene are prominent examples. Beyond fundamental 
interests, these systems might also become instrumental 
in possible future technological advances such as quan- 
tum computing [1]. 

Adding to the excitement is also the possibility to 
simulate a rich variety of models using optical lattices, 
where the dimensionality and interaction parameters are 
highly controllable. With the newly developed method 
of trapping cold dipolar atoms or molecules in such lat- 
tices [2], the particles are subject not only to tunneling 
and contact interaction but also to a long-range interac- 
tion. These systems have been studied intensely (see the 
review [3] and references therein) , and they display a very 
rich phase diagram; Mott or charge-density-wave (CDW) 
states and superfluid and supersolid phases may all be re- 
alized by tuning the strength of the on-site interaction, 
the tunneling amplitude, and the chemical potential. 

In this paper we focus our interest on ultracold bosons 
(atoms or molecules) in a ID optical lattice, considering 
an extended Bose-Hubbard model. The analysis is per- 
formed for a long-range repulsive dipolar interaction, as 
well as for a short-range (on-site plus nearest-neighbor) 
interaction. Various types of CDW phases have been 
investigated in the literature earlier; see, e.g., Refs. [4- 
7]. In this work we add a sequence of exotic members 
to this list by considering a region in parameter space 
in which several CDW's of different symmetries become 



nearly degenerate. While the parameter regions we con- 
sider have been investigated in a few earlier (numeri- 
cal) studies [4, 5], we show that important aspects of 
the problem have been overlooked, in particular, regard- 
ing the entirely new and intriguing excitation structure 
that emerges in this regime. Using a suggestive analogy 
to the physics of non-Abelian fractional quantum Hall 
(FQH) states and an analytical strong-coupling approach 
assisted by numerical calculations, we find a far more 
complex and intriguing phase diagram. First, we consider 
points in parameter space where the problem is exactly 
solvable, leading to a staircase structure of nontrivially 
degenerate gapped CDW ground states. Using intuition 
gathered from recent studies [8, 9] of the Read-Rezayi 
FQH states [10], we find a description of the low-energy 
excitations in terms of domain walls separating degener- 
ate, translationally inequivalent, CDW's. These domain 
walls can be created in inequivalent ways and form a rep- 
resentation of SU (2)fc near lattice filling v = k/2, in close 
analogy with the excitations of the FQH states [9, 10]. 
Within these insulating regimes, fractional charges (or 
fractional particle numbers, in cases where the particles 
are neutral) may be created and manipulated even in the 
system ground state by locally changing the chemical po- 
tential. After having presented this idealized picture, we 
turn to numerical and perturbative calculations and dis- 
cuss what happens away from the exactly solvable points, 
examining different interaction parameters and including 
tunneling. We also point to some experimental implica- 
tions. 

The outline of the paper is as follows. In the next 
section we introduce the effective Hamiltonian describ- 
ing dipolar bosons in ID. The exact ground state of 
the system with neither tunneling nor particle-particle 
interaction beyond nearest neighbors are presented in 
Sec. HI. The low-energy excitations are characterized and 
the analogy of these to the FQH states is outlined. It is. 
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FIG. 1: (Color online) (a) How the CDW states, here for v = fc/2 = 2, can be arranged in a so-called Bratteli diagram. Going 
from left to right, an elementary domain wall excitation is created in each step. A configuration corresponding to the dashed 
red line is shown in (b) and a configuration corresponding to the purple dotted line is shown in (c). In (b) and (c) we show 
only quasi-hole-like domain walls, but quasi-particle-like excitations can be represented analogously. 



furthermore, shown how our results modify when infinite- 
range interaction is included. The proceeding Sec. IV 
considers the effects arising from weak tunnehng between 
the sites. This analysis is performed using perturbation 
arguments as well as numerical calculations, and we find 
that the tunneling in general partially splits the degener- 
acy between the ground states. Finally we conclude with 
a discussion in Sec. V. 



II. EXTENDED BOSE-HUBBARD MODEL 



Vj" > 0, one finds a perfect devil's staircase as a function 
of the chemical potential — this solution is of relevance 
both to the fractional quantum Hall effect [14] and to 
optical lattices [7]. In this paper, we study on-site soft- 
ening to move away from the above convex case and find 
a zoo of new phases and exotic domain walls that should 
be observable experimentally. We will, first, consider a 
truncated Hamiltonian, keeping only the on-site and the 
nearest-neighbor interactions, but we will also show that 
including longer-range terms into the analysis does not 
change our general findings. 



The system of ultracold bosons in a ID optical lattice 
is, after imposing the single-band and tight-binding ap- 
proximations, described by an extended Bose-Hubbard 
Hamiltonian; 

iJ = -t f Qt+i + h.c.) - [L^ fij 

i i j>0 

where t is the amplitude for tunneling between neigh- 
boring sites, n is the chemical potential, U is the on-site 
interaction, and Vj is the energy cost for a pair of parti- 
cles j lattice sites apart. If Vj = Vxjf' is imposed, this 
accounts for a dipole-dipole interaction between the par- 
ticles. The bosonic operator a| (Si) creates (annihilates) 
a particle on site i, while ni = ajoi counts the number 
of particles on site i. In this work we neglect any dipole- 
dipole couplings between neighboring ID tubes. One of 
the two transverse directions can be decoupled by prop- 
erly aligning the dipoles, while coupling in the second 
transverse direction can be minimized by employing a 
sharp harmonic trap [3]. 

As mentioned, the great advantage of utilizing opti- 
cal lattices is the possibility to tune the many system 
parameters [11]; t, U, Vi, and /i can all be adjusted us- 
ing different techniques, like control of the lattice am- 
plitudes, employing Feshbach resonances or alignment of 
the dipoles. It turns out that the sign of the second 
derivative of the interaction, V^' = Vj-i + ~2V} (in- 
cluding VI' = U + V2 — 2Vi), is of crucial importance for 
the low-energy physics of the system [12, 13]. When all 



III. EXACT SOLUTION 

Consider first the Hamiltonian in Eq. (1) with t = 
Vjyi = 0. Assuming that the number of particles is 
allowed to vary, the filling fraction v will be fixed by 
the chemical potential 11. For each such /u, the low- 
energy physics is completely determined by the relation 
between U and Vi or, more specifically, by the sign of 
VI' = U -2Vi. At the point where V(' = 0, i.e., where 
U = 21^1, the following holds: For {k - 1)U < fi< kU, 
k e the filling fraction is v = k/2 and all states with 
periodicity two minimize the energy. In other words, 
the phase diagram is a staircase where the plateau at 
V ^ k/2 contains the CDW states with unit cells [l,k — l], 
I — 0, l,...,fc. Note that there is a (fc + l)-fold de- 
generacy of the ground state. As an example, when 
2>U < n < All, V = 2 and the ground states are the 
degenerate lattice states |04) = [0404...), |13) = [1313...), 
|22) = [2222...), |31) = [3131...), and |40) = [4040...) 
(here, hi\...ni...) = n^\...n,,...)). 

Curiously, precisely these states naturally appear [8, 9, 
15, 16] in the context of the so-called Read-Rezayi [10] 
FQH states. Although the Hamiltonian giving rise to 
the non-Abelian FQH states is starkly different from the 
one studied here, we find that the lowest-lying excita- 
tions of the systems are intimately related. The physical 
reason for this lies within an emergent "generalized ex- 
clusion principle"; the v = k/2 Read-Rezayi states can 
be defined as the unique states within a Landau level 
that vanish as any fc + 1 particles approach the same co- 
ordinate. Analogously, the ground states obtained here 



are the unique states that never have more than k par- 
ticles on two consecutive sites. In fact, when the FQH 
system is put on a thin torus, these two exclusion rules 
become identical (cf. Ref. [9]). Moreover, in both sys- 
tems the lowest-lying excitations consist of fractionally 
charged domain walls between (quasi-)degenerate CDW's 
with unit cells [l,k — I] and [lzLl,k — 1^1], respectively, 
as shown in Figs. 1 and 2 [9]. The fractional charge of 
the excitations is ±|, supposing the particles all carry 
the charge q [17]; to see this, create two identical do- 
main walls by inserting a string of |? ± 1, fc — Z =F 1) in a 
\l,k — I) state, and note that this state has one particle 
more or less than the ground state [9, 18]. For fc > 1 it 
is possible to create domain wall excitations with a given 
charge and position in inequivalent ways, as is demon- 
strated in the left panel of Fig. 1, showing a so-called 
Bratteli diagram [19]. Here, transitions between differ- 
ent CDW states is illustrated by drawing lines connecting 
states with / differing by 1. There are several ways to get 
from the starting to the ending state (the example in the 
figure being |22)), on the way creating a certain number 
of domain wall excitations. Hence, there is a nontrivial 
degeneracy of these. A pedagogical account for the pre- 
cise relation between Bratteli diagrams and topological 
FQH states can be found in Ref. [9] . 

The very same nontrivial degeneracy of the excitations 
is at the heart of the non-Abelian exchange statistics 
in the corresponding FQH states and provides the cru- 
cial non-local degrees of freedom that make these states 
prime candidates in the quest for topologically protected 
quantum computation [1]. However, the analogy between 
the present setup and the two-dimensional FQH states 
has the important limitation that we do not have any 
fractional exchange statistics in our system. Neverthe- 
less, proposals for quantum information processing in ID 
quantum wires have been presented [20]. While we have 
not worked out some similar idea in the current system, 
we may note that the present setup has the nice property 
that more than a single internal electronic level of the 
atom or molecule can be considered. Coupling between 
electronic states then allows for transfer of atoms between 
different atomic lattices, which, in turn, can be spatially 
separated via adiabatic motion. In this respect, single- 
site atoms can be "lifted out" from the lattice [21, 22]. 
Utilizing this possibility might be one way to achieve ex- 
change statistics. 

In order to have a robust realization of the above in- 
troduced fractional domain walls, it is of great interest to 
know whether these types of states can form in a ground- 
state configuration as well. The answer is yes. If the 
experiment is set up so the filling fraction can be forced 
to deviate slightly from half-filling, the fractional charges 
will appear in the ground state. If the filling fraction on 
the other hand is adjusted automatically by the chemi- 
cal potential and a bath of particles, e.g., in the presence 
of a confining harmonic trap, the domain walls can in- 
stead be created by locally varying the chemical poten- 
tial. Single-site control of the chemical potential indeed 
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FIG. 2: (Color online) Schematic demonstration how a local 
change in the chemical potential can give rise to fractional do- 
main walls, causing extended changes in the lattice configura- 
tion. By increasing or decreasing the local chemical potential 
on two neighboring sites, one can create domain walls [(b)- 
(d)] with charge ±g/2, starting from the unperturbed state 
(a). The domain walls shown in this example are closely re- 
lated [8] to the Majorana fermion excitations of the celebrated 
Moore- Read FQH state [23]. 

can be achieved via tightly focused lasers, as was recently 
shown experimentally [21]. This newly developed tech- 
nique for locally controlling /i is most likely superior to 
earlier proposals relying on impurity atoms [24] . We find 
that increasing fi slightly on two neighboring sites will 
tend to create a domain wall with increased particle num- 
ber, while decreasing ^ on two neighboring sites has the 
opposite effect. This is illustrated in Fig. 2. 

In a laboratory setting there is inevitably an external 
trapping potential. This leads to the formation of in- 
compressible domains of different particle densities (each 
at ^' = fc/2 for various fc), separated by relatively sharp 
boundaries. In each of these incompressible domains the 
excitation structure follows the discussion above, thus, 
we expect a given sample to harbor excitations that lo- 
cally mimic those of different FQH states. 

The formation of domain walls is dependent on the 
degeneracy of the different CDW's. A relevant question 
to pose is whether one can tune the interaction so the 
full degeneracy persists also for the more realistic long- 
range dipolar interaction Vj — Vi/j^ (we still treat t — 
0). To answer this question, we note that the energy 
difference between any pair of these states is proportional 
to X]j>i ^2""-!' meaning that if this sum is zero they are 
all degenerate. Inserting the expression for Vj above, one 
finds the requirement 

J>1 J>1 ■' 

= t/- yi^C(3) « t/- 21/1 X 0.9015, (2) 

where C{^) is the Riemann function. Hence, the rela- 
tion between U and Vi, which for the truncated Hamil- 
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tonian was U = 2Vi, has only been slightly shifted, 
and the staircase structure as well as the nature of the 
fractionally charged domain walls remain qualitatively 
the same. However, the domain walls now interact 
weakly due to the long-distance tail of the interaction, 
and the width and positions of the insulating plateaus 
are modified. For the dipolar interaction, the plateau 
corresponding to filling fraction u = k/2 extends over 

uCi^l- Tm^)) < M < c/(ir - 5 + im)^ ""o^s^iy' 

C/(1.1667fc-0.9561) < ii< [/(1.1667fc-0.0439). Between 
these plateaus, which notably do not cover the entire 
range in /i, we argue that the fractional domain walls de- 
scribed above form a devil's staircase of their own — these 
minimally charged excitations repel each other and form 
ordered ground states at any rational filling fraction, just 
like the underlying particles do for a convex interaction 
[7, 13, 14]. 



IV. AWAY FROM THE EXACTLY SOLVABLE 
POINTS 

A. Numerical and perturbative results 

The analysis presented in the previous section is exact 
for f = and when the relative strengths of U and Vi 
are such that the CDW states are exactly degenerate. 
What happens as we relax one of these requirements by 
altering the relationship between U and Vi? We note that 
the domain walls do not interact at the degeneracy point 
U = 2Vi. Slightly away from this point, one finds that 
they interact via a weak linear potential, whose strength 
and sign depends on the energy splittings between the 
CDW's. However, the nature of the domain walls and 
their fractional charge do not change. 

In actual experiments, there is always some finite pos- 
sibility for tunneling between different sites in the lat- 
tice. In this section we present analytical and numerical 
calculations indicating what happens as we include the 
nearest-neighbor tunneling term [t > 0) in the Hamilto- 
nian. Our results are summarized in the phase diagram 
in Fig. 3, which is based on a numerical analysis using the 
Gutzwiller-ansatz method [25] , together with a perturba- 
tive strong-coupling expansion up to third order in t/U . 
Even though it is known that the Gutzwiller mean-field 
approach is not giving a quantitatively correct ground 
state in ID [26], we believe that, together with the per- 
turbative results, it indeed gives a qualitatively accurate 
picture of the phase diagram. 

The degeneracy of the CDW states with unit cells 
[^, k ^ I] will, in general, be split by the introduction 
of the tunneling parameter t. In a strong-coupling ex- 
pansion [26], we include the energy corrections obtained 
using perturbation theory up to third order, and discard- 
ing proportionality factors and redundant constants, the 



energies of the CDW's as a function of I become 

E^l + Ef - [P - kl){U - 2Vi) + 

,9^2 I l(k-l + l) (fc-0(i + i) 1 

\ U(2l-k-l)+Vi(2k-'a + l) ^ U (k-2l-l) + Vi(Al-2k+l) ( ' 

(3) 

assuming Vj>i = 0. This is correct up to third order 
in t/U since all odd-order corrections are zero. Due to 
the form of Eq. (3), it is not possible to make all CDW 
states exactly degenerate simultaneously. However, for 
any fixed value of the tunneling parameter t, one can 
always tune the on-site interaction so CDW states with 
unit cells [l,k — I] and [^ ± 1, A; — / =F 1] become pairwise 
degenerate, and, thus, the fractionally charged domain 
wall excitations remain. To be more specific, since all the 
CDW states are still approximately degenerate for small 
values of t, the low-energy excitations will include all 
types of domain walls contained in the Brattcli diagram; 
see Fig. 1. 

Alternatively, keeping the relationship between U and 
Vi fixed, one can switch between different CDW states 
by varying i, as shown in Fig. 3. Hopping favors states 
with the particles evenly spread out, which means that 
one needs to choose U < 2V\ in order to see the (first- 
order) phase transitions between all CDW states; hence 
Fig. 3 is constructed for U = 1.99Vi. Within the insulat- 
ing lobes in the phase diagram, the energy difference be- 
tween any two CDW states at the specific filling fraction 
is independent of jj,. Consequently, the lines separating 
two adjacent CDW phases are vertical. The positions of 
the vertical lines in the figure have been calculated using 
Eq. (3). 

The remaining structure of the phase diagram has been 
obtained numerically using the Gutzwiller-ansatz wave 
function 

i \ji=0 / 

where \n)i is the n-atom Fock state at site i, fn^ is the 
corresponding amplitude, and A'^ is a truncation number 
[N — 30 in our calculations). We restrict the numerics 
to the most relevant situations consisting of solutions of 
periodicity two, and, accordingly, the lattice occupations 
fn^ may be separated into two sublattices A and B [4]. 

The values of /l'^'^^ are obtained by minimizing the en- 
ergy functional E = {'$\H\'^) self-consistently. For the 
given parameters we find a set of different supersolid (SS) 
phases [27] surrounding the insulating CDW's. These 
SS's are characterized by nonzero 0q = (na) for the two 
lattices a = A, B, together with a nonzero order param- 
eter A(/) — (pA — 4'B- The different CDW's as well as the 
SS's are both separated by first-order phase transitions. 
In the insulator regime it follows that the transitions be- 
tween CDW crystalline phases have to be of first order 
since must be discontinuous at the transition. We 
believe that the same mechanism drives the transitions 
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between SS phases, and, hence, the transition is, again, 
of first order. In particular, the SS's can be seen as par- 
tially melted CDW phases, and the different SS phases 
support distinct underlying CDW states. Further away 
from the insulating regimes, the SS phases turn into pure 
superfluid phases {4'a ~ <Pb 7^ 0) via second-order phase 
transitions. In order to clearly distinguish between the 
different SS phases. Fig. 3 is constructed based on calcu- 
lations of the derivative dA(j)/dt of the order parameter. 




t/U 

FIG. 3: (Color online.) Phase diagram for U = 1.99Vi, 
Vj>i = 0, obtained by employing a Gutzwiller ansatz, to- 
gether with perturbation theory up to third order in t/U. The 
yellow vertical lines represent the phase transitions between 
the different CDW states, and the diagram also displays the 
transitions between the insulating CDW regions, the super- 
solid (SS) phases and the superfluid (SF) region. The col- 
orscale marks only the different phases. 



B. Experimental implications 

The different regimes in the phase diagram, presented 
in Fig. 3, should be accessible experimentally (see discus- 
sion in Ref . [4] ) , opening up the possibility of seeing new 
types of CDW phases in a ID lattice. Another prospect 
discussed in our study concerns the possibility of vary- 
ing the chemical potential locally to create fractionally 
charged domain walls. A suggestion might be to switch 
on the disturbance in /i adiabatically, using the technique 
of Ref. [21], and then suddenly turn it off again to fol- 
low the system evolution as the domain wall delocalizes, 
giving a measure of its effective mass or the so-called 
Lieb- Robinson bound [28]. Similar measurements were 
indeed recently performed in a system following a quan- 
tum quench [30]. As for this experiment, the dynamical 
delocalization can, in principle, be tracked down with 
the help of individual-site fluorescence methods, recently 
implemented in various experiments [29]. One impor- 
tant observation is that the above on-site measurements 



are limited to detect the "parity," meaning the even or 
odd number of atoms. Even though there are sugges- 
tions regarding how to go beyond this restriction [31], we 
note that whenever a domain wall separates two different 
CDW phases, the parity differs between the two phases. 
Thus, even if the measurement is capable of detecting 
only the parity, this should already contain information 
of the propagation of the excitations. An alternative way 
to create such a nonequilibrium state could be to drive 
the system nonadiabatically through one of the phase 
transitions separating the various CDW's and thereby to 
create defect excitations. 



V. CONCLUDING REMARKS 

We have studied ultracold bosons in a ID optical lat- 
tice, assuming an extended Bose-Hubbard model. For 
particular choices of the strength of the on-site interac- 
tion U, we find a nontrivial increased ground-state degen- 
eracy, both for a truncated Hamiltonian and for a long- 
range dipolar interaction between the particles. In the 
weak-tunneling regime, the ground states are different 
CDW states with periodicity two. Within the degenerate 
regime it is possible to create fractionally charged domain 
walls between the degenerate vacua by varying the chem- 
ical potential locally in the lattice. Including tunneling 
in the analysis, we presented a phase diagram achieved 
using perturbation theory together with a numerical sim- 
ulation assuming a Gutzwiller ansatz, and the type of 
phase transitions between the different CDW states, as 
well as between supersolid and superfluid regimes, were 
identified. We also pointed out that all experimental 
techniques necessary to verify our results already exist. 

In this work we also presented an intriguing relation- 
ship between ultracold bosons in ID and the 2D FQH 
effect. Wc showed that on-site softening leads to an 
emergent "generalized exclusion principle" that forbids 
fc -|- 1 particles to occupy two consecutive sites, in strik- 
ing similarity to the clustering conditions on the non- 
Abelian Read Rezayi FQH states that dictates how the 
wave-function amplitude must vanish as fc -I- 1 approach 
the same coordinate. In fact, in both systems this prop- 
erty can be taken as the defining property of the ground 
states (as well as the quasihole excitations). We also 
showed that the structure of low-energy excitations, built 
from "fractional domain walls," is identical to that of 
the Read Rezayi states. Moreover, we submit that this 
"generalized exclusion principle" should be much easier 
to directly observe in the ID dipolar systems discussed 
here than in the context of non-Abclian FQH states 
where such concepts are frequently used; see, e.g., Refs. 
[8, 9, 15, 16]. 

We also want to stress that the analogy between ultra- 
cold bosons in ID and the 2D FQH effect also holds in 
the case of simpler, Abelian FQH states; the solution of 
the very same problem giving rise to the devil's staircase 
for dipolar bosons (without on-site softening) [7] was ear- 
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lier used to provide a microscopic understanding of the 
self-similar phase diagram observed in the lowest Landau 
level [14]. 

We hope that this analogy can spur further cross- 
fertilization between the two fields and we note that it 
is possible to obtain a whole plethora of nontrivially de- 
generate CDW ground states by making other choices for 
various V^' [12]. 
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